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Sum rules for twist-3 distributions are reexamined. Integral relations between 
twist-3 and twist-2 parton distributions suggest the possibility for a 5-function at 
x = 0. We confirm and clarify this result by constructing h^ and h^ (quark-gluon 
interaction dependent part of h^) explicitly from their moments for a one-loop 
dressed massive quark. The physics of these results is illustrated by calculating 
h]^(x,Q 2 ) using light-front time-ordered pQCD to 0(ctg) on a quark target. 

1. Introduction 

Ongoing experiments with polarized beams and/or targets conducted at RHIC, 
HERMES and COMPASS etc are providing us with important information 
on the spin distribution carried by quarks and gluons in the nucleon. They 
are also enabling us to extract information on the higher twist distributions 
which represent the effect of quark-gluon correlations. In particular, the twist- 
3 distributions gxiXjQ 2 ) and 1il(x,Q 2 ) are unique in that they appear as 
a leading contribution in some spin asymmetries: For example, gx can be 
measured in the transversely polarized lepton-nucleon deep inelastic scattering 
and Iil appears in the longitudinal-transverse spin asymmetry in the polarized 
nucleon- nucleon Drell-Yan process 1 . The purpose of this paper is to reexamine 
the validity of the sum rules for these twist-3 distributions. 

A complete list of twist-3 quark distributions is given by the light-cone 
corelation functions in a hadron with momentum P, spin 5* and mass M: 



= 2 [ 9l (x, Q 2 )p»(S ■ n) + g T (x, Q 2 )S^ + M 2 g 3 (x, Q 2 )n»{S ■ »)] , (1) 
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/ ^e iX *(PS\ma^i~f 5 i;(\n)\ Q 2\PS) = 2 [h^x, Q 2 )(S^ - S^)/M 
+h L (x, Q 2 )M{p»n v - p v n^){S ■ n) + h 3 (x, Q 2 )M(S^ - , (2) 



dX 
— e 
2tt 



lXx (PS${0)iP{\n)\ Q 2\PS) = 2Me{x,Q 2 ). (3) 



The light-like vectors p and n are introduced by the relation p 2 = n 2 = 0, 
n+ =p- = 0, P» =p^ + ML n n an d S» = (S ■ n)p^ + (S ■ p)n» + S£. The vari- 
able x G [— 1, 1] represents the parton's light-cone momentum fraction. Anti- 
quark distributions gi t T, 3 (x, Q 2 ), h\^^{x, Q 2 ) are obtained by replacing ip by 
its charge conjugation field C^> T in (l)-(3) and are related to the quark distri- 
butions as gi,T,3(x, Q 2 ) = gi,T,z{-x, Q 2 ) and h ljTt3 (x, Q 2 ) = -h liTi3 (-x, Q 2 ). 
The sum rules in our interest are obtained by taking the first moment of the 
above relations. For example, from (1), one obtains 

(ps\mri5m\Q*\ps) w 

= ^J dx [ 9l (x,Q 2 )pV(S ■ n) + g T (x,Q 2 )S*+ M 2 g 3 (x , Q 2 )n» (S ■ n)] . 

From rotational invariance, it follows that the left hand side of (4) is propor- 
tional to the spin vector and thus gi,r, 3 {x, Q 2 ) must satisfy 

J dx gi (x,Q 2 ) = J dxg T {x,Q 2 ), (5) 

J dx 9l (x,Q 2 ) = 2 J dxg 3 (x,Q 2 ). (6) 
The same argument for (2) leads to the sum rule relations for hi^, 3 {x, Q 2 ): 
J dxh 1 (x,Q 2 )=J dxh L {x,Q 2 ), (7) 

J dxh x {x,Q 2 ) = 2 J dxh 3 (x,Q 2 ). (8) 

The sum rule (5) is known as Burkhardt-Cottingham sum rule 2 and (7) was 
first derived in Refs. 3 ' 4 . Since the twist-4 distributions g 3 , h 3 are unlikely to be 
measured experimentally, the sum rules involving those functions (6) and (8) 
are practically useless and will not be addressed in the subsequent discussions. 
As is clear from the above derivation, these sum rules are mere consequences of 
rotational invariance and there is no doubt in their validity in a mathematical 
sense. However, if one tries to confirm those sum rules by experiment, great 
care is required to perform the integral including x = 0. In DIS, x is identified 
as the Bjorkcn variable xb — Q 2 /2P ■ q and x — corresponds to P ■ q — > oo 
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which can never be achieved in a rigorous sense. Accordingly, if tiL(x,Q 2 ) 
has a contribution proportional to 5(x) and h\(x,Q 2 ) does not, experimental 
measurement would claim a violation of the sum rule. 

In this paper we reexamine the sum rule involving the first moment of the 
twist-3 distribution Iil{x,Q 2 ). In particular we argue that Iil{x,Q 2 ) has a 
(^-singularity at x = 0. Starting from the general QCD-based decomposition 
of Hl, we show that it contains a function /i™ which has a <5(x)-singularity. 
In Sec. 3, we construct Hl for a massive quark from the moments of h\ at the 
one-loop level and show that h\ also has an (5(a;)-singularity, which together 
with the singularity in /i™ gives rise to a S(x) singularity in itself. In Sec. 4, 
we perform an explicit light-cone calculation of in the one-loop level to 
confirm the result of the previous sections. Details can be found in Ref. 5 . 



2. (5(:c)-functions in hz(x, Q 2 ) 

The OPE analysis of the correlation function (2) allows us to decompose 
hL{x,Q 2 ) into the contribution expressed in terms of twist-2 distributions 
and the rest which we call h\{x,Q 2 ). Since the scale dependence of each 
distribution is inessential in the following discussion, we shall omit it in this 
section for simplicity. Introducing the notation for the moments on [—1,1], 
M n [hL\ = J_j dx x n h,L(x), this decomposition is given in terms of the moment 
relation 1 : 



M M = — 
M [h L ] = M [/ii 



M n [hi] + 



n m q 
71 + 2M" 



M„_i[5i] +M n [h 3 L ], 



(n>l) (1) 
(2) 



with the conditions 



M [h 3 L } = M 1 [h 3 L } = 0. 



(3) 



By inverting the moment relation, one finds 



h L (x) = hY w (x) + h?(x)+hl(x) 



= < 



2x 



J 

J x 



dy 



m„ 



-2x 



y 2 M 

hi(y) 



2x 



Jx y A 



+ hl(x) 



dy 



V 



M 



2x 



9i(y) 



+ h 3 L (x), 



(4) 

(x > 0) 
(5) 
(x < 0) 
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where the first and second terms in Eq. (4) denote the corresponding terms in 
(5). In this notation the sum rule (2) and the condition (3) implies a 

M o [/#]=0. (6) 

Naively integrating (5) over x for x > 0, while dropping all surface terms 6 one 
arrives at dxh^x) = dxhi(x) + dxh\(x) and likewise for dxhL{x). 

Together with (3), this yields J^dxh^ = j\dxhi. However, this procedure 
may be wrong due to the very singular behavior of the functions involved near 
x = 0. Investigating this issue will be the main purpose of this paper. 

We first address the potential singularity at x = in the integral expression 
for h™(x) in (5). In order to regulate the region near x = 0, we first multiply 
h™(x) by x 13 , integrate from to 1 and let /3 — > 0. This yields 

/ dxV£(x) = ^ lim P { dxx^gUy) = ^<?i(0+), (7) 

while multiplying Eq. (5) by \x\P and integration from —1 to yields 

j\ dxh?(x) = lim/3 |° dx\x\^ 9l (y) = -^i(O-), (8) 

where we have assumed that <?i(0±) is finite. Adding these results we have 

|° dxh?{x) + jf dxhJ£{x) = £± ( 9l (0+) 9l (0-)) . (9) 

Eq. (6) and the fact that, in general, linx^o 9i{%) — g\ (—x) ^ 0, b imply 

TTt 

hl£{x) = hT(x) reg - ^ ( 3l (0+) - ffl (0-)) S(x), (10) 

where h™(x) reg is defined by the integral in (5) at x > and x < and is 
regular at x = 0. Eq. (10) indicates that Hl has a S(x) term unless h\(x) has 
a (5(x) term and it cancels the above singularity in h™(x). 

Eq. (10) demonstrates that the functions constituting Hl{x) are more sin- 
gular near x = than previously assumed and great care needs to be taken 
when replacing integrals over nonzero values of x by integrals that involve the 
origin. In particular, if h^x) itself contains a 8{x) term, then (2) implies 

[ (h L (x)-h 1 (x))+ [ (h L (x) - h^x)) ^0, (11) 

a More precisely, the original OPE tells us Mo[h A L + h™] = 0. But as long as gi(0±) is finite, 

which wc will assume, this is equivalent to stronger relations (3) and (6). 

b For example, dressing a quark at 0(as) yields gi(0+) ^ and gi(0— ) = <?i(0+) = 0. 



5 



and, since h\(x) is singularity free at x = 0: 

/O- rl rl 

dxhi(x) + / dxhi(x) = / dxh\(x). (12) 
-l 7o+ 7-1 

Accordingly an attempt to verify the "/ii,-sum rule" 3 would obviously fail. 

However, in order to see whether the 5(x) identified in (10) eventually 
survives in 1il(x), we have to investigate the behavior of h? L (x) at x = 0. To 
this end we will explicitly construct hi(x) for a massive quark to 0(as). 



3. Hl(x, Q 2 ) from the moment relations 

In this section we will construct hi(x,Q 2 ) for a massive quark (mass m q ) to 
0(as) from the one-loop calculation of M„[/i 3 ]. 



h L (x,Q z ) = h%'(x) + -^C F ln^h^(x), (13) 

All 777 , . 



where the scale Q is introduced as an ultraviolet cutoff and the Cp = 4/3 is the 
color factor. h V ^ W,3,m ^ ,1 \x) are defined similarly. gf\x) = h^\x) = 5(1 — x) 
gives hf\x) = 5(1 — x), as it should. One loop calculations for g\(x) and h\(x) 
for a quark yield the well known splitting functions 7 ' 8 : 

h^(x) = ^^- + \5(l-x). (15) 

Inserting these equations into the defining equation in (5), one obtains 

h™ W{1) (x) = Zx + Ax\n-^, (16) 

h ? {1) (x) = 2 - 4a;lni^ - 3 + 3x - i<f(a;). (17) 

In the first line of (17), the term (3x — §(5(1 — x)) comes from the self-energy 
correction, i.e. from expanding M = m q 1 + f^CF§m^? m Eq. (4) ; an d 

— |#(a;) = — ^gi(0+)5(x) in h™^\x) accounts for the second term on the r.h.s. 
of (10). We also note that hY W ^ does not have any singularity at x = and 
satisfies f} dx \i^ W ^ (x) = f} dx (x) as it should. 

h^\x) can be constructed if we know the purely twist-3 part h^ 1 \x) at 
the one-loop level. One- loop renormalization of hi was completed in 9 and the 
mixing matrix for the local operators contributing to the moments of hi(x, Q 2 ) 
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was presented. We obtain for the moment of the quark distribution 5 



J dxx n h 3 L {1 \x) = 



6 



n + 1 



1 

2 + 2 



(18) 



3( 1) 

for n > 2. From this result, together with (3), we can construct /i L (x) as 



h$ 1 \x)=3-6x+±6{l-x)-±6{x) 



(19) 



We emphasize that the —l/2S(x) in (19) is necessary to reproduce the n = 
moment of h^ 1] {x). From (16), (17) and (19), one obtains 



h L (x,Q 2 ) = S(l-x) + ^ln^C F 



• (20) 



We remark that the above calculation indicates that the 8{x) term appears not 
only in h™ but also in h\. Furthermore they do not cancel but add up to give 
rise to — 5(x) in Hl(x,Q 2 ) itself. In the next section we will confirm Eq. (20) 
through a direct calculation of hi,(x, Q 2 ) for a quark. 



4. Light-cone calculation of fiL(x,Q 2 ) 

In order to illustrate the physical origin of the 5(x) terms in Hl(x) and to de- 
velop a more convenient procedure for calculating such terms, we now evaluate 
Iil{x) using time-ordered light-front (LF) perturbation theory. The method 
has been outlined in Ref. 10 and we will restrict ourselves here to the essen- 
tial steps only. There are two equivalent ways to perform time-ordered LF 
perturbation theory: one can either work with the LF Hamiltonian for QCD 
and perform old-fashioned perturbation theory 10 , or one can start from Feyn- 
man perturbation theory and integrate over the LF-energy fc~ first. In the 
following, we will use the latter approach for the one- loop calculation of Hl(x). 

In LF gauge, A + = 0, parton distributions can be expressed in terms 
of LF momentum densities (fc + -densities). Therefore, one finds for a parton 
distribution, characterized by the Dirac matrix T at 0(as) and for < x < 1 

Mx Mp) ru {p) = -i^v) J ^^k^lk^Te^ 

xS^x-^PjD^ip-k), (21) 
where D^(q) = g^ u - Mbd^nSa j s the gauge field propagator in LF 



gauge, and is a light-like vector such that nA = A + ~ (A + A 3 ) j\[2 for 
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any four vector A^. The k integrals in expressions like Eq. (21) are performed 
using Cauchy's theorem, yielding for < k + < p + 

dk~ 1 1 



/ 



2tt ( fc 2_ TO 2 + ie ) 2 (p-k) 2 +ie 

1 1 1 k,-oo 1 1 



(22) 



(2fc+) 2 2(p+ - k+) ( m'+kj ( P ^-k^ \ 2 2p+ kl ' 

^ 2fc+ 2(p+-fc+) ^ 

where we used k + = xp + . In order to integrate all terms in Eq. (21) over k~ , 
Cauchy's theorem is used to replace any factors of k~ in the numerator of Eq. 
(21) containing k~ by their on-shell value at the pole of the gluon propagator 

k~ — jfe- = v- - (p± ~ k±)2 (23) 

In the following we will focus on the UV divergent contributions to the parton 
distribution only. This helps to keep the necessary algebra at a reasonable 
level. We find for < x < 1 to 0(a s ) 

M^> = f>l»|-^, (24) 

where the usual +-prescription for [ 1 _ 1 a; ] + applies at x = 1, i.e. [iz^jr. = jzj 
for x < 1 and J* dx [ 1 _ 1 rE ] + = 0. Furthermore, hi,{x) = for x < 0, since 
anti-quarks do not occur in the O(as) dressing of a quark. In addition to Eq. 
(24), there is also an explicit 5(x— 1) contribution at x = 1. These are familiar 
from twist-2 distributions, where they reflect the fact that the probability to 
find the quark as a bare quark is less than one due to the dressing with gluons. 
For higher-twist distributions, the wave function renormalization contributes is 
t^Cf In Qztj8(x—1). The same wave function renormalization also contributes 
at twist-3. However, for all higher twist distributions there is an additional 
source for 5(x — 1) terms which has, in parton language, more the appearance 
of a vertex correction, but which arises in fact from the gauge-piece of self- 
energies connected to the vertex by an 'instantaneous fermion propagator' 
t^jt. For gr{x,Q 2 ) these have been calculated in Ref. 10 where they give an 

additional contribution — f^Cp In %?$( x ~ 1)> i- e - the total contribution at 
x = 1 for gr(x, Q 2 ) was found to be ff-Cp In ^i(5(x— 1). We found the same 
S(x — 1) terms also for Hl(x, Q 2 ). c Combining the S(x — 1) piece with Eq. 



c In LF gauge, different components of the fermion field aquire different wave function renor- 
malization. However, since all twist-3 parton distributions involve one LC-good and one 
LC-bad component, one finds the same wave function renormalization for all three twist-3 
distributions. 
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(24) we thus find for < x < 1 

h L (x, Q 2 ) = 5(x - 1) + ^C F In 9" 



[1 - x]+ 2 

Comparing this result with the well known result for h\ ' 



+ h(x-l) 



hl (x,Q 2 )=6(x-l) + ^\n^C F 
2tt nig 



o*(*-l) 



[1 - x}+ 2 



(25) 



(26) 



one realizes that 



lim / * dx [h L {x, Q 2 ) - h x {x, Q 2 ] = ^ In ^rC F ? 0, (27) 

i.e. if one excludes the possibly problematic region x = 0, then the fei-sum 
rule 3 is violated already for a quark dressed with gluons at order 0(as)- 

In the above calculation, we carefully avoided the point x = 0. For most 
values of fc + , the denominator in (21) contains three powers of k~ when fc~ — > 
oo. However, when fc+ = 0, fc 2 — ra 2 becomes independent of k~ and the 
denominator in (21) contains only one power of k~ . Therefore, for those terms 
in the numerator which are linear in k~ , the ^"-integral diverges linearly. 
Although this happens only for a point of measure zero (namely at k + = 0), 
a linear divergence is indicative of a singularity of h L (x,Q 2 ) at that point. d 
To investigate the k + w singularity in these terms further, we consider 

/(fc+k ± )= [dk-— k -L— - (28) 

7 (A; 2 — to 2 + ie) 2 (p — k) A + is 

k~ + (k~ - fc 



^ (fc2_ TO 2 +ie )2[(p_ fc )2 +ie ] -/ca n .(fc + ,ki) + / sm .(fc+,k ± ), 

where the 'canonical' piece f ca n. is obtained by substituting for fc~ its on 
energy-shell value fc~ = p~ — ^2(p+-k+) t tnc vame at tne P°l e at (p — fc) 2 = , 
Eq. (23)]. For fc+ = 7^ 0, it is only this canonical piece which contributes. 
To see this, we note that fc~ — fc~ = — 2 (p+-\+) > anc ^ therefore 



kj.) = J dk 



1 



(fc 2 — m 2 + ie) 2 (p — fc) 2 + ze 
1 f 1 

= 2(p+ - fc+) 7 ^ (fc 2 - to 2 + ie) 2 ' (29) 

Obviously 12 J (2fc+fc--k^ -mg+ie) 2 = for fc + 7^ because then one can 

m 2 _[_J^ 2 J£ 

always avoid enclosing the pole at k~ = — q 2k + by closing the contour in 



d Note that the divergence at fc+ = p+ is only logarithmic. 
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the appropriate half-plane of the complex k — plane. However, on the other 
hand jd 2 k L [kl _^J: ml+ . e)2 = and therefore 

fsUk\^) = ^^l. (30) 



2p+ k 2 ± 



m 
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Upon collecting all terms oc k in the numerator of (21), and applying (30) to 
those terms we find for the terms in Iil(x, Q 2 ) that are singular at x = 

h L ,Hnix, Q 2 ) = ~ ^ %C F 5{x). (31) 
z7r ra q 

Together with Eq. (25), this gives our final result for Hl, up to 0(as), valid 
also for x = 

a s n , Q 2 
— C F In — 

2tt mi 



-S(x) + ^ + 1^-1) 



.(32) 



h L {x, Q z ) = 5(x - 1) + ^ C F In ^ 

~i q 

As expected, ft./, from Eq. (32) does now satisfy the /il-suhi rule, provided of 
course the origin is included in the integration. 

This result is important for several reasons. First of all it confirms our 
result for tiL(x,Q 2 ) as determined from the moment relations. Secondly, it 
provides us with a method for calculating these S(x) terms and thus enabling 
us to address the issue of validity of the naive sum rules more systematically. 
And finally, it shows that there is a close relationship between these S(x) terms 
and the infamous zero- modes in LF field theory 13 . 

Ref. 11 , where canonical Hamiltonian light-cone perturbation theory is used 
to calculate Hl(x) and, for x ^ the result obtained in Ref. 11 agrees with 
ours which provides an independent check of the formalism and the algebra. 
However, the canonical light-cone perturbation theory used in Ref. 11 is not 
adequate for studying the point x = 0. From the smooth behaviour of /i_l(x) 
near x — the authors of Ref. 11 conclude that the sum rule for the parton 
distribution Hl(x) is violated to 0(as). Our explicit calculation for Hl{x) not 
only proves that the sum rule for hi,{x) is not violated to this order if the 
point x = is properly included, but also shows that it is incorrect to draw 
conclusions from smooth behaviout near x — about the behaviour at x = 0. 



5. Summary 

We have investigated the twist-3 distribution Hl(x), and found that the sum- 
rule for its lowest moment is violated if the point x = is not properly included. 
For a massive quark, to 0(as) we found 



h L (x, Q 2 ) = S(x 



l) + ^ln^ 
m 2 



2tt 



-S(x) 



[1 - x}. 



\S(X 



1) 



(1) 



10 



At O(as), h,L,(x, Q 2 ) does not satisfy its sum rule if one excludes the origin from 
the region of integration (which normally happens in experimental attempts to 
verify a sum rule). Of course, QCD is a strongly interacting theory and parton 
distribution functions in QCD are nonpcrturbative observables. Nevertheless, 
if one can show that a sum rule fails already in perturbation theory, then this is 
usually a very strong indication that the sum rule also fails nonperturbatively 
(while the converse is often not the case!). 

From the QCD equations of motion, we were able to show nonperturbatively 
that the difference between Iil(x,Q 2 ) and h\{x,Q 2 ) contains a 5{x) term 

[h L (x, Q 2 ) - hl(x, Q 2 )] singular = (. gi (0+, Q 2 ) 91 (0-, Q 2 )) 5(x). (2) 

Since gi(0+,Q 2 ) — 51 (0— ,Q 2 ) = lim^o 9i( x , Q 2 ) — gi(x.Q 2 ) seems to be 
nonzero (it may even diverge 6 ), one can thus conclude that either hi,{x,Q 2 ) 
or h\{x,Q 2 ) or both do contain such a singular term. 

We checked the validity of this relation to O(as) and found that, to this 
order, both h\ and contain a term oc S(x). We also verified that even 
though the sum rule for /il(x) is violated if x = is not included, it is still 
satisfied to O(as) if the contribution from x = (the 8{x) term) is included. 
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